Adiabatic Geometric Phase for a General Quantum State 
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A geometric phase is found for a general quantum state that undergoes adiabatic evolution. For 
the case of eigenstates, it reduces to the original Berry's phase. Such a phase is applicable in both 
linear and nonlinear quantum systems. Furthermore, this new phase is related to Hannay's angles 
as we find that these angles, a classical concept, can arise naturally in quantum systems. The results 
are demonstrated with a two-level model. 
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Consider a quantum system that depends on some ex- 
ternal parameters R. We are interested in the evolu- 
tion of its quantum state when the parameters R change 
slowly along a closed path (see FigQ. For an eigen- 
state, such an adiabatic evolution accumulates a geomet- 
ric phase, which reflects the system geometry with the 
parameter space R. It is now known as Berry's phased). 
Although the same geometry is also embedded in the 
adiabatic evolution of a general quantum state, includ- 
ing non-eigenstate, it is not clear how a geometric phase 
can be defined for a general quantum state to extract the 
system geometry. 

In this Letter we introduce a geometric phase for a gen- 
eral quantum state. In the case of eigenstates, this phase 
can be reduced naturally to the original Berry's phase. 
From this reduction, one gains some fresh perspective on 
the familiar concept of Berry's phase. 

The geometric nature of this new phase is fully ex- 
plored. Similar to Berry's phase, it can be regarded as 
the geometric part of a total phase associated with the 
adiabatic evolution. In another perspective, it can also 
be regarded as a part of an Aharonov-Anandan (AA) 
phasejj, 0- Moreover, the new phase is found to be 
related to Hannay's angles^ in a derivative form. Al- 
though Hannay's angles have long been considered as the 
counterpart of Berry's phases in classical mechanics, we 
find that they can be defined naturally in quantum sys- 
tems, and be calculated from the new phase. In this 
sense, this relation unifies two very different concepts, 
Berry's phase and Hannay's angles. 

Interestingly, this new geometric phase is also applica- 
ble in nonlinear quantum systems. Unlike in linear quan- 
tum systems where one may understand the adiabatic 
evolution of a non-eigenstates in terms of eigenstates 0, 
it is impossible to do the same in nonlinear quantum 
systems due to the lack of the superposition principle. 
Therefore, the introduction of this new phase provides 
a unique and powerful tool to study the adiabatic evo- 
lution of a general state in nonlinear quantum systems. 
Bose-Einstein condensates of dilute atomic gases are an 



excellent example of nonlinear quantum systems 0, Q- 




FIG. 1: Adiabatic evolution of a quantum state along a closed 
path in the parameter space of R. In general, the quantum 
state does not come back at the end of the evolution, that is, 
the difference between the start state ip and the end state ip' 
is more than a phase. The disks represent the dynamical evo- 
lution of a quantum state at a given R; their different shapes 
symbolically illustrate the change of the system Hamiltonian 
with R. The areas of these disks remain the same, reflecting 
that the actions i" are conserved. 

In the following discussions, we will proceed with non- 
linear quantum systems since linear quantum systems 
can be regarded as their special cases. 

We consider an TV-level quantum system governed by 
a general nonlinear Schrodinger equation (H = 1), 



i-\i,) = H{r j ^ j ;Rm 



(i) 



Here \ip) — (ipi, ip2, ' ' ' > V'jv) is the wave function with ipj 
being its jth component over an orthonormal basis; the 
vector R represents all the system parameters subject to 
adiabatic change. We assume that the system is gauge 
invariant since it is the case of most physical interest. 
When H is independent of ip* and ipj, Eq. (JJJ is the 
usual linear Schrodinger equation. 
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It is well known that the quantum system governed by 
Eq. Q mathematically has a canonical classical Hamilto- 
nian structure (e.g., see Refs.@,[||). That is, one can find 
a Hamiltonian H(ip*, ipj] R) such that Eq. Q is a set of 
equations of motion when 7i serves the classical Hamil- 
tonian with Poisson brackets {ipj,ipk} = ifijk- When 
the Hamiltonian is bilinear, Ti = {ip\H\ip), the system is 
linear. With this Hamiltonian structure, the system of 
Eq.([Q) can be classified into integrable or non-integrable 
in the classical sense. All linear quantum systems are 
integrable classically. 

We focus on the case that Eq.Q is integrable. In this 
case, the system at a given R has N constants of motion, 
J = {7i, I2, ■ • • , In}- They are called actions, whose con- 
jugate variables are angles — {9i,9 2 , ■ ■ ■ , 9jy}, which 
change with time linearly as 9 = 0q + ujt, where ujj ~ 
dTL/dlj. Therefore, at a given R, the wave function can 
be expressed as a function of {/, 6}, \ip) = \ifj(I,0; R)). 
Such a parameterization of the wave function in terms 
of {/, 6} is also very convenient for slowly changing 
R. In such an adiabatic evolution, the actions I are 
still conserved [Toj while the angles change as 0(t) = 
Oq + J u>(R)dt + a(t), where a are Hannay's angles|4j. 

With the above observations, we are now mathemat- 
ically ready to introduce a geometric phase for the adi- 
abatic evolution of a general quantum state |i/>}, where 
the parameters R change slowly along a closed path O, 
as shown in Fig^ One feature immediately stands out: 
Except eigenstates, the quantum state generally does 
not come back to the original state at the end of evo- 
lution even though the system Hamiltonian recovers its 
original form. Moreover, the difference between the ini- 
tial state and the ending state depends on the choice of 
the initial state. These features make it difficult to de- 
fine Berry's phase and its various generalizations for this 
evolution 0, 

To circumvent the obstacle, we use an averaging tech- 
nique as in the definition of Hannay's angles[4|. For a 
quantum states — \ip(I,0; R)) adiabatically evolving 
along the path O, we introduce a geometric phase 7 as 

l[(D) = £ dRWI > 9; m ik m ' 6; R)) ' (2) 

where the integration over 6 is to average over all pos- 
sible quantum states with the same actions 7 at a given 
R. We emphasize that the integration over R in Eq.© 
is done for a fixed 6. This new geometric phase 7 is the 
same for all the quantum states that has the same ac- 
tions I; therefore, it may also be regarded as a geometric 
phase for an adiabatically evolving manifold. To illus- 
trate and further explore this new phase 7, we start with 
some special cases and end with a demonstration with a 
nonlinear two-level model. 

For a quantum system starting at an eigenstate 



\E n (R)) defined by 

H(R)\E n (R)) = E n (R)\E n (R)), (3) 

it evolves dynamically as \ip(t)) — e~^^\E n (R)) when 
R changes adiabatically. In this special case, there is 
only one action, the norm I = (E n (R)\E n (R)) , whose 
corresponding angle is the phase 9 = (3. Plugging this 
\%p) into Eq.© and noticing that the partial derivative 
over R does not act on 9 and /, we obtain the phase 7 
for eigenstate \E n ), 

7n (0) = £dR{E n (R)\i-^\E n (R)) , (4) 

which recovers the original Berry's phase [lj. Note that 
Eq. is valid in both linear and nonlinear quantum sys- 
tems, indicating that the original definition of Berry's 
phase can be directly borrowed for nonlinear quantum 
systems if only eigenstates are considered. 

We turn to a general quantum state in linear quantum 
systems. In this linear case, we can expand the evolving 
quantum state in terms of the eigenstates, 

JV 

|^(t)> = ^a n (t)| J B n (Jl)>. (5) 

71=1 

According to the quantum adiabatic theorem^3]j the oc- 
cupation probabilities of different eigenstates \a n \ 2 are 
adiabatic constants. In fact, they are actions /„ = \a n \ 2 
when the system is regarded mathematically as a classical 
Hamiltonian system; their corresponding angle variables 
9 n are the phases of a n 's. With these in mind, com- 
putation of Eq.® with the state (J5J is straightforward. 
We find that the off-diagonal terms are zeros after the 
averaging, and the geometric phase 7 is 

JV 

7(0H5>„| 2 7„(0), (6) 

n=l 

where 7„ is the Berry's phase of eigenstate \E n ). There- 
fore, in linear quantum systems, the phase 7 is just 
a weighted summation of the Berry's phases of all the 
eigenstates involved. Interestingly, this kind of weighted 
summation of Berry's phases has already been applied in 
calculating transverse fore on a quantized vortex ll 31 and 
the anomalous Hall conductivity of ferromagnets 14] . 

We have so far illustrated the geometric phase 7 for 
some simple examples and shown clearly how the phase 7 
as defined in Eq. |j2l are related to the well-known Berry's 
phase. In the following, we are going to examine the new 
phase 7 in a general setting and derive for it a different 
expression. These efforts reveal that the phase 7 can be 
regarded as a geometric part of the total phase of an 
adiabatic evolution, similar to Berry's phase. 

Imagine an adiabatic evolution of a general quantum 
state \ip) along a close path O. Without worrying about 
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whether it is geometric or not, we can always introduce 
a phase for such an evolution, 



(3 = l J dt (^\j t \^) = i f o ^\ d \^) ' ( 7 ) 

where T is the total evolution time. We call this phase 
the total phase of the adiabatic evolution of \ip). We 
expand the integrand as 

(</>(!, 9; R)\dMI, 0; R)) = <V>I^> • dl + 



+ {^\ — \i,).de + 



l dR 



\ijj)-dR. (8) 



For an adiabatic evolution, the actions are conserved so 
the term involving I is zero. Plugging it back into Eq.Q) 
and averaging it over the initial angles 6 , we obtain 



& = 



d0 

(2ir) N 



= I ■ j dO + 7 (G), (9) 

where we have introduced a new notation (3 to stand 
for the averaged total phase f3. We may also call j3 the 
total phase of an adiabatically evolving manifold. In the 
derivation, we have noticed for 7 that the integral over 
is the same as over a at any given R. We have also 
used that 



(10) 



is the action of the motion associated with angle 8j . Note 
that this expression is also the Aharonov-Anandan (AA) 
phase of the cyclic state, describing the motion associated 
with the angle variable 9j . 

The equation 10 shows that the averaged total phase 
(3 has two parts: a dynamical part involving action-angle 
variables and a geometric part that is exactly our new 
phase 7. For the special case of eigenstate = \E n ), 
since there is only one non-zero action / = (E n \E n ), we 
have 

p = p={E n \E n ) [ u> n dt + ln {0). (11) 
Jo 

It recovers the well-known fact that Berry's phase is the 
geometric part of the total phase of an adiabatic evolu- 
tion of an eigenstate 1] . 

We re- write Eq.© and obtain another expression for 
the geometric phase 7, 



dO. 



(12) 



This complicated expression turns out to be easier to be 
implemented numerically; we will use it to compute the 
results shown in Fig [21 

One can repeat the derivation from Eq.(0 to Ea. l)12|l 
in the projective Hilbert space. In this situation, the 
total phase Eq.Q becomes an AA phase, and the phase 
7 can then be regarded as the geometric part of the AA 
phase associated with the parameter space. 

There is another interesting angle looking into the new 
phase 7. As already mentioned, the system described in 
Eq.Q has a canonical classical Hamiltonian structure. 
With this classical structure, we can introduce naturally 
Hannay's angles, a classical concept, into quantum sys- 
tems. Following the expression for Hannay's angles in 
Ref. 0, 0, 0] , we find that these angles in our quantum 
system Q are related to the phase 7 by 



a(0) = -Jj7(O). 



(13) 



In linear quantum systems, these Hannay's angles dif- 
fer from Berry's phases of eigenstates only by a sign, 
a n = — 7„, according to Eq.@). Note the relation ltT^|) 
is very similar in form to the semiclassical relation be- 
tween Hannay's angles and Berry's phases derived in 
Ref.0 [13 although they are two very different rela- 
tions. So far, it is not clear why the similarity |is|. 

Finally, we demonstrate the geometric phase 7 with a 
nonlinear two-level model as given by 



. dipi 
1 dt 
. d(f 2 
1 dt 



c\<pi\ 



</?2 



X-iY 
2 

X + iY 



■ 



-fi ■ 



(14) 
(15) 



This simple model can be used to describe the Josephson 
effect of Bose-Einstein condensates residing in a double- 
well potential|3, Il9j| . The complex coupling constant, 
as denoted by X and Y, can be realized in experiment 
through phase imprinting on one of the two wells |20j. 

We first look at the geometric phase 7 for eigenstates 
of this nonlinear quantum system. This is to find all the 
eigenstates \E n (R)) for a closed path O, and use them to 
calculate the phase with Eq. Q . Both steps could be done 
numerically; fortunately for this simple case, analytical 
results can be obtained. When the path is restricted on 
the unit sphere X 2 + Y 2 + Z 2 — 1, we find that the phases 
for these nonlinear eigenstates are 



1(0) 



rf(R + cqz) -dS 
9S=0 (c V + Z) 2 (c^ + Z) 



(16) 



where z is the unit vector along the z-axis and r] is one 
of the real roots of 



cV + 2cZrf + (1 - c 2 )r] 2 - 2cZr) - Z 



0. 



(17) 



Different real roots r\ correspond to different eigen- 
states. It is clear that Eq. l|17|l can have more than two 
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are the analytical results of Eq. (|16fl . 
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FIG. 2: Geometrical phases 7 of the nonlinear two- level model 
1141151 : (a) the lower eigenstate and (b) the non-eigenstate 
near the lower eigenstate. The non-eigenstate has an action 
(or A A phase) of I = 0.005. The insert shows the closed path 
O used which is a circle on a unit sphere with Z fixed. The 
circle is traversed with a rate of 0.001. The results for eigen- 
states are compared to the analytical expression of Eq. ll6H 
(denoted by diamonds): a very good match is found. 



real roots, indicating that there can be more than two 
eigenstates[2lj. Here we limit ourselves to the situations 
where Ea. 1)17(1 has only two real roots. For the path O 
that is a circle with a fixed Z , the geometric phase in 
Ea. (|16fl becomes 7 = (1 — 77)77. The diamonds in Fig|3 
are calculated with Eci. (|16J) . showing how 7 for the lower 
eigenstate changes with Z. 

For a general quantum state, we have to resort to nu- 
merical means. The path O is picked to be a circle with 
fixed Z. We then solve Eas. l|14ll5|l numerically after 
choosing a changing rate v = 0.001 for the parameters 
R = {X, Y, Z}. The evolving states are recorded and 
used to compute the phase 7 with Ea. (fl"2")l . where the av- 
eraging is done for different initial states with the same 
action (or AA phase) /. Results for c = 0.05 are plotted 
in Fig|2 showing how the phase changes with Z. Com- 
putation is done for both eigenstate and non-eigenstate 
and the results (solid lines) are compared to the phases 
for the linear case c = 0.0 (dashed lines). The changing 
rate of R (v = 0.001) is slow enough to be considered as 
adiabatic. This is witnessed by the good fit between the 
solid line and the diamonds in Fig[3fa) as the diamonds 



In summary, we have introduced a geometric phase 
for a general adiabatically evolving quantum state. The 
new phase to certain extent unifies two different concepts, 
Berry's phase and Hannay's angles. It is very interest- 
ing to find potential applications for this new geometric 
phase while its properties are being further explored. 
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